Density expectation value of two independent interacting Bose-Einstein condensates 



in 

O 

o 

(N 
00 



5h 
^3 



c3 



o 

(N 
> 

m 

m 

o 
in 
o 

-I— > 



i 

c 

o 
o 



X 
S3 



Hongwei Xiong 1,2,3 , Shujuan Liu 1,2,3 , Mingsheng Zhan 1,2 
1 State Key Laboratory of Magnetic Resonance and Atomic and Molecular Physics, 
Wuhan Institute of Physics and Mathematics, Chinese Academy of Sciences, Wuhan 430071, P. R. China 
2 Center for Cold Atom Physics, Chinese Academy of Sciences, Wuhan 430071, P. R. China and 
3 Graduate School of the Chinese Academy of Sciences, P. R. China 
(Dated: February 2, 2008) 

After removing the double- well potential trapping two initially independent Bose condensates, the 
density expectation value is calculated when both the exchange symmetry of identical bosons and 
interatomic interaction are considered. After the overlapping, it is shown that there is a nonzero 
interference term in the density expectation value. This nonzero interference term physically arises 
from the exchange symmetry of identical bosons and interatomic interaction which make two initially 
independent condensates become coherent after the overlapping. It is found that the calculated 
density expectation value with this model agrees with the interference pattern observed in the 
experiment by Andrews et al (Science 275, 637 (1997)). 
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The coherence property plays an essential role in the 
wave nature of Bose-Einstein condensates (BECs), and 
it has been investigated intensively [lj after the exper- 
imental realizations of BECs in dilute gases. For two 
coherently separated Bose condensates, it is not surpris- 
ing that there is a clear interference pattern when the two 
condensates are allowed to overlap. In the celebrated ex- 
periment by Andrews et al Q, however, high-contrast 
fringes were observed even for two completely indepen- 
dent condensates at an initial time. This experimental 
result shows clearly that there is a coherence property 
after the overlapping of two independent interacting con- 
densates. 

In many literature (see for example 0,3 and references 
therein), upon expansion, it is shown that there is no in- 
terference term in the density expectation value for two 
initially independent condensates, and thus the observed 
high-contrast fringes [2j for two independent condensates 
were interpreted |3| with the aid of the high-order corre- 
lation function P (r, r', t) (which is an oscillation function 
of r — r') and quantum measurement process. The con- 
tradiction between the vanishing interference term in the 
density expectation value and the observed high-contrast 
fringes was also discussed based on several theories such 
as the stochastic simulations of the photon detection for 
atoms [3 , the expansion of Fock state by the linear super- 
position of coherent states , and the continuous mea- 
surement theory ■ In addition, the nonlinear effects in 
the interference pattern of two coherently separated Bose 
condensates were investigated in Refs ja, . 

In the present work, we calculate carefully the density 
expectation value for two initially independent conden- 
sates by including the interatomic interaction. It is found 
that for the case of two initially independent condensates, 
upon expansion, there is a nonzero interference term in 
the density expectation value when interatomic interac- 
tion and the exchange symmetry of identical bosons are 
taken into account carefully! After removing the double- 
well potential trapping the two initially independent con- 
densates, we give the theoretical result of the density ex- 



pectation value which agrees with the interference pat- 
tern observed in . 

For the condensate trapped in a magnetic trap, it can 
be divided into two condensates by a blue-detuned laser 
beam. The combined double-well potential takes the 
form 
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where the second term represents the external potential 
due to the laser beam. 

After the evaporative cooling for Bose gases trapped 
in the double-well potential, if the intensity of the blue- 
detuned laser beam is sufficiently low so that the tun- 
neling effect is obvious, the two condensates can be re- 
garded to be coherently separated. For two coherently 
separated condensates, every atom is described by the 
following normalization wavefunction 



4>c (r, t) 



W^i (r, t) 



/VN, (2) 



where <f>\ (r, t) and (f>2 (r, t) are the normalization wave 
functions accounting for the two condensates. N — 
Ni + N2 is the total number of particles. Before the 
overlapping of the two condensates, the average particle 
numbers in each condensate are Ni and N2, respectively. 
After removing the double- well potential, the density ex- 
pectation value is then 

n c (r,t) = N\a c \<f> 1 {r,t)\ 2 + 2b c xRe(<f)* 1 (r,t)(P 2 (r,t)) 



+c c \<p2 (r,t)r 



(3) 



where a c = Ni/N, b c = y /N 1 N 2 /N and c c = N 2 /N. 
The second term in the above equation accounts for the 
interference effect when there is an overlapping between 
two condensates upon expansion. 

If the intensity of the blue-detuned laser beam is suffi- 
ciently high so that the tunneling effect can be omitted, 
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the two condensates can be regarded to be completely 
independent. In this case, the number of particles Ni 
and N 2 in each of the two condensates is fixed £|. To 
investigate clearly the role of the exchange symmetry of 
identical particles and interatomic interaction, we calcu- 
late the density expectation value rid (r, t) directly from 
the many-body wave function. Based on the well-known 
exchange symmetry of identical bosons [lOj, the many- 
body wavefunction of the whole system is given by 



^n 1 n 2 Oi,r 2 , 
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iVi!iV 2 ! 
NI 



^2P[<pi (ri,i)-->i (r Nl ,t)x 



(r Nl+ i,t) • • - 02 (r Nl+N2 ,t)] , 
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where P denotes the Nl/(NilN2\) permutations for 
the particles in different single-particle state (f>i or 
<f>2- The above expression of the many-body wave- 
function is based on the quantum mechanical prin- 
ciple 10]. A n is a normalization factor to assure 
Xl*WiJV 2 {ri,r 2 , - ■ ■ ^N,t)\ 2 dr 1 dr 2 - ■ - dr N = 1. A n is 
determined by the following equation: 



A ,, 



where 
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C(t)= / <t>i(v>t)4>* 2 (T,t)dV=\C(t)\e i 



(5) 
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In this Letter, to give a concise expression for various 
coefficients such as A n , we have introduced the rule 0° = 
1. 

From the above many-body wavefunction, the exact 
density expectation value takes the following form: 

n d (r,t) = N J^* NlN2 (r,r 2 ,---,r N ,t) x 

V Nl N 2 (r,r 2 ,---,rjv,f) dr 2 ■ --dr N 

1 • 2b d x Re (e**" 1 ' 



= NB, 
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where B n = A^Ni\N 2 \/N\, and the coefficients in the 
above equation are given by 
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For two independent ideal condensates, before the 
overlapping of the two condensates, £ (t = 0) = 0. Based 
on the Schrodinger equation, it is easy to verify that af- 
ter the double-well potential separating the condensates 
is removed, one has £ (t) = at any further time. Thus 
bd = 0, and the density expectation value is given by 



n d (r, t) = Ni |0! (r, t)\ 2 + N 2 |0 2 (r, t)Y 



(11) 



In this situation, the interference term is zero in the den- 
sity expectation value. 

However, if the interatomic interaction is considered, 
£ (t) can be a nonzero value, and this will lead to a quite 
different result compared with the ideal condensates. As- 
suming that g = 4irh 2 a s /m with a s being the scattering 
length, after removing the double-well potential, the en- 
ergy of the whole system is given by 
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i<3 



The second term in the above equation represents the 
interaction energy of the whole system by using the well- 
known two-body pseudopotentials. 

By using the ordinary action principle and the energy 
of the whole system, one can get the following coupled 
evolution equations for tfii and <f> 2 : 



ih- 



ih- 



dt 



dt 



1 SE 

1 SE 
N~ 2 W 2 



(13) 
(14) 



where 8E/8(f>l and SE/S4> 2 are functional derivatives. 
With the above coupled equations, one can understand 
that C (t) becomes nonzero after the overlapping be- 
tween two initially independent condensates for g being 
a nonzero value |ll| . 

Although generally speaking, |C(^)| is much smaller 
than 1 because (j>i(/) 2 is an oscillation function about the 
space coordinate, nevertheless, a nonzero value of £ (t) 
will give significant contribution to the density expecta- 
tion value for large Ni and N 2 . In Fig. 1(a), we give the 
relation between bd/cd and £ for Ni = N 2 = 10 3 based 
on the numerical calculations of Eqs. (JSJ, 10 and 110|l . 
The relation between bd/cd and Ni — N 2 for £ = 0.001 
is shown in Fig. 1(b). It is shown clearly that even for 
ICI being the order of N± , the value of bd can not be 
omitted, and thus the interference term in Eq. (JJJ will 
play an important role in the density expectation value. 
Generally speaking, for JVi |C| >> 1 and 7V 2 |C| >> 1, 
one has b d / \fa d c d « 1. 

Generally speaking, increasing the particle number will 
enhance the effect of the interference term in the density 
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FIG. 1: Based on the numerical calculations of Eqs. (SJ, © 
and HUH , shown in Fig. 1(a) is the relation between bd/cd 
and while shown in Fig. 1(b) is the relation between bd/cd 
and JVi = N2 for ( = 0.001. It is shown clearly that the 
interference term can play an important role in the density 
expectation value for \(\ being larger than JV-T . 



FIG. 2: Shown is the evolution of the density expectation 
value calculated using the experimental parameters in Ref. 
0. It is shown that there is a clear interference pattern after 
the overlapping of two independent interacting condensates. 
Shown in the inset is the evolution of |C(^)| for the experi- 
mental parameters in Ref. |2j. 



expectation value. Based on Eqs. (| 1 31) and l|14l) . increas- 
ing the coupling constant g will have the effect of in- 
creasing £(i). Together with the relation between bd/cd 
and £ illustrated in Fig. 1(a), this shows that increasing 
the interatomic interaction will increase the effect of the 
interference term. 

In the experiment of Ref. |2j, N = 5 x 10 6 condensed 
sodium atoms were confined in a magnetic trap with 
lu x = 2ix x 18 Hz, uj y = uj z = 2ir x 320 Hz. A blue- 
detuned laser beam of wavelength 514 nm was focused 
into a light sheet with a cross section of 12 /im by 67 
/im. The long axis of the laser beam was perpendicular 
to the long x— axis of the condensate. For a laser power 
of 14 mW, the barrier height is about 1.4 /iK, which 
is much larger than the chemical potential /i = 0.03 fiK. 
For this laser beam, the two condensates can be regarded 
to be independent because they are well separated and 
the tunneling effect can be omitted (See also 0). With 
these experimental parameters and the s— wave scatter- 
ing length a s = 2.75 nm, the initial profile of the two 
condensates is shown in Fig. 2. After the double- well 
potential is removed, the evolution of the density expec- 
tation value rid-x (x, t) = J rid (r, t) dydz (in unit of N/2) 
is given in Fig. 2 through the numerical calculations of 
Eqs. 0, {T2Jl and {Tl}. We see that there is a clear in- 
terference pattern in the density expectation value which 
agrees with the experimental result. Shown in the inset 
is the evolution of |C| for these parameters. We have also 
verified in the numerical calculation that |£| is always 
zero for two ideal condensates. For the expansion time 
of 40 ms, the numerical result of \(\ shows that bd/cd~ 1. 

Based on the evolution of the density expectation value 
shown in Fig. 2, we see that the overall width (about 



300 fim) of the ultracold gases in x— direction does not 
increase obviously. This is due to the fact that the initial 
density distribution is cigar-shaped, and thus the expan- 
sion in £— direction is very slow, while the expansion in y 
and z directions is much quick. In the experiment of Ref. 
|l2j |. one can see clearly that there is no obvious expan- 
sion in the long a;— axis for cigar-shaped condensate. For 
two initially independent condensates, when the double- 
well potential is switched off, one should note that in the 
regime close to x = 250 /im shown in Fig. 2, the ul- 
tracold gases expand rapidly in a;— direction because in 
this regime the ultracold gases have higher kinetic en- 
ergy. Thus, although the total width of the system in 
x— direction does not increase obviously, the rapid ex- 
pansion in the central regime will lead to the overlapping 
between two initially independent condensates, and re- 
sults in the interference effect. The overall width of the 
system shown in Fig. 2 is smaller than the experimental 
result of about 500 fxm 0. This difference may come 
from the expansion of thermal cloud in this experiment 
0. After 40 ms expansion, the numerical result in Q for 
two coherently separated condensates also showed that 
the overall width of the system is about 300 /xm. 

The essential reason for the emergence of the interfer- 
ence term of two initially independent condensates lies 
in that because of the exchange symmetry of identical 
bosons and interatomic interaction, the two initially in- 
dependent condensates will become coherent after the 
overlapping between the two condensates. The physi- 
cal mechanism of this interaction-induced coherence can 
be also understood based on the second quantization 
method. For two initially independent condensates, the 
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quantum state of the whole system can be written as |4j 
\N 1 ,N 2 ) = 7 ^_(a\r(ai) N >\0), (15) 

where C n is a normalization constant to assure 
(Ni, N 2 \Ni, N 2 ) = 1. a\ (aj) is a creation operator which 
creates a particle described by the single-particle state </>i 
(</>2) in the left (right) condensate. The density expecta- 
tion value is then 

n d (r,t) = (N 1 ,N 2 ,t\^\N 1 ,N 2 ,t). (16) 

Here ^ (x,i) is the field operator. The operators <Si 
and a 2 can be written as ai = J^tfildV and a,2 = 
J ^(fi^dV , respectively. By using the commutation re- 
lations of the field operators (x, £) , \t (y, t)] = and 

(x, i) , vfrf (y, £)] = <5 (x — y). It is easy to get the com- 
mutation relation [aijcrf,] = C*- We see that ai andaJ, are 
not commutative any more for J tfii^dV being a nonzero 
value. The non-commutative property between ai andaJ, 
plays an essential role in the emergence of the interfer- 
ence term. Our calculations show that the result given by 
Eq. is the same as the result given by Eq. which 
is obtained based on the many-body wavefunction. We 
have also proven that the evolution equations for </>i and 
4>2 based on the many-body wavefunction are the same 
as the results based on the second quantization method 

El. 

One can get a further physical picture for the emer- 
gence of interference pattern of two initially independent 
condensates through a general investigation for the case 
of N x |C| » 1, N 2 |CI >> 1 and N x ~ N 2 . Based on 
our numerical calculations, these conditions are satisfied 
for the case of 40 ms expansion in the experiment of Ref. 
0- When these conditions are satisfied, NB n ad ~ Ni, 
NB n b d « y / N 1 N 2 and NB n c d w N 2 . Thus, the density 
expectation value can be approximated very well as 

n d (r,i)~$*(r,i)$ e (r,i), (17) 



where the effective order parameter $ e (r, t) is 

$ e (r, t) = (r, i) + v^Vbe*^ (r, t) . (18) 

In particular, after removing the double-well potential, 
the evolution of the effective order parameter $ e (r, t) can 
be described very well by the following Gross-Pitaevskii 
equation El 



iH—^~-—V 2 <S> e +g\<S> e \ 2 <S> e . (19) 



We see that the emergence of the effective order pa- 
rameter gives us strong evidence that the coherence is 
formed between two initially independent condensates, 
and thus results in the emergence of high-contrast inter- 
ference fringes. 

In summary, upon expansion, we calculate the density 
expectation value of two initially independent conden- 
sates. Quite different from the current beliefs, it is found 
that there is a nonzero interference term when the inter- 
atomic interaction and the exchange symmetry of iden- 
tical bosons are both considered carefully. In fact, it is 
well-known that the interaction plays an essential role in 
the formation of the order parameter of Bose-condensed 
gases, i.e. the formation of a stable coherent property. 
Here, we provide an example in which the interaction in- 
duces coherent evolution between two initially indepen- 
dent condensates. To investigate more clearly the inter- 
ference pattern due to the measurement process alone, 
an experimental investigation of two independent ideal 
condensates would be very interesting, because there is 
no interference term for ideal condensates in the den- 
sity expectation value. In the last few years, the rapid 
experimental advances of Feshbach resonance where the 
scattering length can be tuned from positive to negative 
make this sort of experiment be feasible. 
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